مبانی ترمودینامیک در پایتون 


(معادلات حالت و مدل‌های ضریب فعالیت) 


۱۱0۲ 5086۲ )۲,( 


6010 5۱۷ 0( 


مترجم وتهیه کننده : علی زین الدینی 


کتاب مبائی ترمودینامیک و پایتون (جلد اول) یکی از کتاب های مورد استفاده برای محققان شیمی و علاقه مندان 
حوزه برنامه نویسی پایتون است. در اف کتاب سعی شده است تا مبانی ترمودینامیک در قالب کدهای بر نامه نویسی 
پایتون در اختیار علاقه مندان قرار بگیرد. در جلد اول درمورد معادلات حالت 9 ضریب فعالیت بحت شده انستت و با 
استفاده از کتابخانه 167100 در پایتون این مباحث به زبان برنامه نویسی جهت محاسبات بهتر کد نویسی گردیده 
تمامبی مطالب استفاده شده ترجمه وب سایت زیر است. همچنین نحوه استفاده از این کتابخانه نیز شرح داده شده 


0 و( 


این کتاب توسط علی زین الدینی دانشجو رشته شیمی فیزیک و محقق در حوزه سیستم های شیمیایی و 
الکتروشیمیایی و توسعه دهنده پایتون تهیه گردیده است. برای ارتباط با مترجم و تهیه کننده میتوانید با ایمیل زیر 
مکاتبه داشته باشید. همچنین لینک صفحه اینستاگرامی ما نیز به آدرس زیر خواهد بود 


ون, قورع ۵ 118۷1996 ملع 112 :۴2211 
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۱ ص جح 4 


نقل قول ی ها ی واه هو ها 
شاخص ها 9 جداول مه و هسام اه مج ده مه وج عم ه قع وم اه دعوم ماه وه هدعم و دودعم با اوه و هه ونم دزیم مه امه هم و هه مدوم وج 
۱) مقدمه ای بر معادلات حالت مکعبی وا تا موب هنعط ماه ام دس مخ وگ مق و وب 

اع) کار نا اشنا علض ۳ 


۲-۱) تعادل اجزای خالص 1 


)کار پاتمتخاوظ ها ی 


۵-۱) تعادل مخلوط ره 
۶-۱) استفاده از واحدها با معادلات حالت مکعبی هب سوه موه یو مهم 


۲) مقدمه ای بر مدل های ضریب فعالیت که همه وی وم واه و6 هه وا هو و وهی 
ت۱۷ 9[ 


۲-۳) نکاتی در مورد عملکرد هب کر وود ی دی کب هه 


۲-۴) ویژگی های دیگر 7 


۵-۳۲) هویت های ضریب فعالیت ی ی و ی هیبعت او دق ی رو 


بط 


و هت هو ی ید هس و ۳ 


۱ 


هه ی کت 


0 ۳ 


0 کتابخانه ملفه ترمودینامیک و تعادل فاز طراحی مهندسی شیمی (51)) 


نحوه نصب کتابخانه 1067۳۳0 در بایتون 
آخرین نسخه 1067۳00 را از لینک زیر دریافت کنید 

00/0( 
اگر پایتون را با پیپ نصب کرده اید. به سادگی آن را با استفاده از: 

ممحصعط الماعصز ونم ٩‏ 

از طرف دیگر اگر از 60802 به عنوان مدیریت بسته خود استفاده می کنید. می توانید به سادگی 1۳67۳00 را از کانال 
002-60 در محیط خود نصب کنید: 

مصصمطا 0۲26]-02صمع ع- ااقافصا قجمع ٩‏ 
برای دریافت نسخه )1 فرمان زیرا اجرا کنید 


نع .مصصعط /11 طحعله) نمی حنطانع//:انع عصمل ]نع ٩‏ 


آخرین کد منبع 
آخرین نسخه توسعه منابع 1۳67۳00 را می توانید در این آدرس دریافت کنید 


1۱/60/19۵5[ ظ حاعله) /ع۵ع .ما تزع //:۲)۵5 


اطلاعات مجوز 

برای اطلاعات در مورد شرایط و ضوابط استفاده از این نرم افزار. و سلب مسئولیت از کلیه ضمانت هاء به 0۳5۳.06[ 
که آدرس آن در زیر قرار دارد مراجعه کنید. اگرچه 1067۳00 نیازی به مجوز آن ندارده اما اگر برای شما راحت است. 
لطفاً در صورت استفاده از ۳6۲00 در کار خود ذکر کنید. همچنین در نظر داشته باشید که هر تغییری که ایجاد می 
کنید برای جامعه و توسعه دهندگان مفید است. 


۱۱۱۱ /ماومهه/۵۱09/مصصعط /۱ 01 ظ ماعله۲) /عر۵ع .ما تزع //:۲)۵5ظ 


گزارش اشکال 
برای گزارش اشکالات. لطفاً از 1۳2616 ع3 دمصع1 در آدرس زیر استفاده کنید: 


۵[۱/)۵0/1690۵5ظ حاعله) /عم .ما تزع //:۲)۵5 


۳ ص دح 4 


اگر سوّال بیشتری در مورد استفاده از کتابخانه دارید. با نویسنده در 2011.07 22160.۸070767.13۴611)0۵) تماس 


۳ 


نقل قول 


برای استناد به 10671۳00" در نشریات از عبارت زیر استفاده کنید: 


061713 62۵1 0015 


(بزراب۲ 


شاخص ها و جداول 


) 61 


۱۵115 )2016-2021( ۰. ۵110 : 


۳۲۱01۳66۳13۲0 ۳۲653 0۴ ۷ 


]0۵8 2۳80 11ظ وه 021 


د 0 06104 ۵۴ ]0۵۲8۵۵8۵ 


۰ 161۳0 611ظله 6۵8/081 . طتاط] 7/03 : 1۳09 


برای دسترسی به شاخص ها و جدول های این کتابخانه می‌توانید به آدرس های مشخص شده مراجعه کنید 


شاشفر شا 
شاخص های ماژول 


۳ ص غ جح ه 


) مقدمه ای بر معادلات حالت مکعبی 


معادلات حالت مکعبی» مدل‌های ترمودینامیکی سازگار و نسبتاً سریعی را برای مواد شیمیایی و مخلوط‌های خالص 
ارائه می‌کنند. آنها معمولاً برای نمایش گازها و مایعات استفاده می شوند. 


فرم سه پارامتری عمومی به شرح زیر است: 


( 6611 ۱:۷۹ 
ع 4 ۵۷ +۲۷۶ 6 ۷ 


این اساس پیاده سازی در ترمو را تشکیل می دهد. دو رابط مجزا ارائه شده است. 067۳00.609] برای مدل‌سازی جزء 
خالص و 067100.609_70] برای مدل‌سازی چند جزئی. اجزای خالص کمی سریعتر از مخلوط های چند جزئی هستند. 
زیرا قوانین اختللاط واندروالس که به طور معمول مورد استفاده قرار می گیرند. عملیات 2 ۳ برای محاسبه (0 
انجام می دهند: 
رز(۵۵) زمرت ظ 8 > :060 
۳ 

جمزاام5_حم_وعز۳۳0261_ام5 است) است که سایر ویژگی‌ها ممکن است به آن‌ها بستگی داشته باشند. و محاسبه 
حجم مولی با توجه به یک جفت ورودی (۳ ,1) «کل یک زیرماژول 0611۳00.605_۷0101706]) این موضوع را مورد بحث 
و اجرا قرار می دهد. هر دوی این محاسبات در زمان ثابت هستند. بنابراین سربار آنها برای اجزای خالص و مخلوط‌های 


۱-۱ کار با اجزای خالص 
ما می‌توانیم از رابط 0205 (مخفف «معادله عمومی حالت مکعب ») با هر جزء یا معادله پیاده‌سازی‌شده حالت 
استفاده کنیم. اما برای سادگیء 9 -هگزان با 05 «مومز-۳۵۵ استفاده می‌شود. دمای بحرانی آن ۵۰۷,۶ کلوین. 


هنگام تنظیم اولیه یک شی 02005 وضعیت باید همراه با ثابت های بحرانی مشخص شود. ما در اینجا از > ۴۰۰ و 
2 166 استفاده می کنیم: 


امد ۵0 ۶0 << 
(1۳6عط ,1400 ,۰.29175 ۳۵9220 ,03025000۰0 ,۰6 ۲62501) ۲58 < ووه <<< 
(0. ۶1000000<ظ ,۲2۶400۰0 ,2975 1۵920ه , ۰0 3025000عظ ,6 ۲62507) ۲58 ووم <<< 


۵ ص و< 4 


رشته _ 760۳ _ طوری طراحی شده است که تمام ورودی های شی را نشان دهد. ما می توانیم راه حل های حجم را 


با ویژگی 127_۷010۳065 بررسی کنیم: 


<<< 6095 . 127 7010۵ 
)0۰0001560731847856, 0۰002141876816741, 0۰00091929 5474982( 


ه 1۳ ۱ .۰ 


صفت فاز در این حالت مقدار "1/2" خواهد داشت. در غیر این صورت پا "ع خواهد بود. 
6 5 <<< 


۱1/9 


ویژگی های اساسی محاسبه شده در مقداردهی اولیه فستفیها خواص هستند و می توان به آنها دسترسی داشت: 


خواص مایع مانند در انتهای نام خود "1_" دارند و "ع_" در انتهای خواص گاز مانند است. 


<<<۵09 , ۲ 060 1 


/ 0 
06۳0 09.5 <<< 
8 6۵- 
1 0۳ 0۴ . وه << 
3. 2889501 


تمام محاسبات در 067100.605] و 67100.605_101] بر اساس مولار است. وزن مولکولی هرگز ارائه یا مورد نیاز نیست. 
همه خروجی ها در واحدهای ٩‏ پایه (6 ۳۵2 103 مول. و غیره) هستند. این توسعه را به طور قابل ملاحظه ای ساده 
کرد. برای کار با واحدهای مبتنی بر انبوه» از رابط فاز استفاده کنید. رابط‌های 6۳0.605 و «067100.605_۳0] قبل 
از رابط فاز توسعه داده شده‌اند و دارای برخی ویژگی‌ها هستند که در اینترفیس فاز نمایش داده نشده‌اند. سایر ویژگی 
ها یا به عنوان روش هایی که به آرگومان نیاز دارند پیاده سازی می شوند. یا ویژگی های پایتون که دقیقاً مانند 
خواص‌ها عمل می کنند اما نتایج را به سرعت محاسبه می کنند. به عنوان مثال. فوگاسیتی فاز مایع 2201071 با 
ضریب انبساط هم‌باریک گاز (فشار ثابت) خواص هستند. 


<<<۵05 ۰۴1920111۷ 1 
2 9 5 

<<< 6095.6] 
0-0 9 


۶ ص جح ه 


تعداد بسیار زیادی از این ویژگی ها وجود داردء زیرا بسیاری از آنها مشتقاتی هستند که شرایط مشابهی دارند. فهرست 
کامل در اسناد ۳05/)) موجود است. تماس‌های کمتری وجود دارد که دما را می‌گیرد. مانند 11720 که گرمای تبخیر 


جسم را در دمای مشخص محاسبه می کند: 


(300) ۲1720 . وهوه <<< 
0(0> 


هنگامی که یک شی ایجاد شد. می توان از آن برای نمونه سازی اشیاء جدید 0/۳05 در شرایط مختلف بدون تعیین 


مجدد ثابت های بحرانی و سایر پارامترهای مورد نیاز استفاده کرد. 


(15عظ ,۰0 ۲<300) 0۵9.10 <<< 

(100000۰0عظ ,۲۶300۰0 ,۰.2975 ۳9۵9220ه , ۰.0 3025000ععظ ,۰.6 ۲0۵2507) ۲5٩‏ 
(15ع<ظ ,۷162 ۵09.10 <<< 

(۰0 ۷00 ,00000۰0 1عظ ,۰.2975 20و۳9۵ه , ۰:0 3025000ععظ ,۰.6 ۲02507) ۲5٩‏ 
( ۲۶300 ,۷162 ۵09.10 <<< 

(۰0 ۷۶00 ,۲۶300 ,۰.2975 ۳۵9220ه ,۰.0 3025000عظ ,۰.6 ۲02507) ۲5٩‏ 


همانطور که در مثال های بالا مشاهده شد. از هر دو 1 ۴ ۷ می توان برای مشخص کردن وضعیت جسم استفاده 


کرد. متغیرهای ورودی شی ذخیره می شوند و می توان آنها را با 51216_50605 بررسی کرد: 


09: 9131 , وم <<ح 
1000000۰01 ۲5۲۶ ,400.0 ۲۲۲۶ ) 


بخش های جداگانه معادله عمومی مکعبی نیز ذخیره می شود. می توانیم از آنها برای بررسی اینکه معادله فشار برآورده 
شده است استفاده کتیم! 


5 1۳۳001۲ 605 . ۵۳0 درگ <<< 

+ 46112 . ومع*1 ۷ ووم + 1*۲2 ۰7 605) /21۳02 ۵ ووه - (ظ, 1-605 ۰7۷ 605) /1 .3۳*605 <<< 
(103 ۵0951 . وم 

1000000 0 

+ 0612 . 605* 6095۰7 + 2** ۰7 605) /21012 ,ومع - (ط. 05ع- ۰۲ 605) ۰.۲ ۲5*605 <<< 
(10۳3 6091 , وه 

1000000 0 


۷ ص جح هه 


توجه داشته باشید که از آنجایی که نقاط شناور کاملا دقیق نیستند. ممکن است برخی خطاهای کوچک نشان داده 


شود. اما دقت زیادی برای به حداقل رساندن ان انجام شده است. 


مقدار ثابت گاز مورد استفاده (1/)01*16 ۱۴۴۶۲۶۱۸۱۵۲۲۴ ۸:۲) است. این نزدیک به دقت کامل اعداد ممیز شناور 
است. اما نه کاملاٌ اکنون این یک مقدار دقیق است که به عنوان "تعریف " در سیستم 81 استفاده می شود. توجه داشته 
باشید که سایر پیاده‌سازی‌های معادلات حالت ممکن است از مقدار کامل ثابت گاز استفاده نکنند اما نویسنده اکیدا 
توصیه می کند هر کسی که پیاده‌سازی ۳05 خود را بنویسد از ثابت گاز کامل استفاده کند. این به نتایج قابل تعویض 


بیشتری اجازه می دهد. 


تعانال الح‌ای خالص 
در ادامه با همان حالت و مثال قبلی دو راه حل از معادله حالت موجود بود. با این حال. مگر اينکه دمای دقیق ۴۰۰ 
کلوین و فشار 1 مگاپاسکال روی خط اشباع باشد. هميشه یک حالت ترمودینامیکی پایدار دیگر وجود دارد. ما باید از 
انرژی آزاد گیبس خروجی استفاده کنیم تا مشخص کنیم کدام حالت پایدارتر است. برای یک جزء خالص, حالتی که 
خروج انرژی آزاد گیبس را به حداقل می رساند» پایدارترین حالت است. 

(1۲6عط , ,1۲400 ,2975 0ح2و۳9۵ه ,۰.0 3025000عظ ,۰.6 507ع۲6۵) ۲58 < ووه <<< 


<<< 095.6 060 1, 605 ۰.06060 5 
)-2872.498434, -973 ۰ 5198207( 


به راحتی می توان دید فاز مایع پایدارتر است. این میانبر استفاده از انرژی آزاد گیبس خروجی فقط برای اجزای خالص 
تا باه کارهات کا ال سح کار آیده ان انشفادهمی کفنه متفر است, عیار کال آیخ ات که کرام از 
انرژی آزاد واقعی گیبس را به حداقل برساند. متد 076_502016_0۳256 این بررسی را انجام می دهد و "1" یا "2" را 
برمی گرداند: 


هط 52016 1۵1۲6 . وم <<< 
بآ 


برای یک جزء خالص. یک خط تعادل بخارمایع درست تا نقطه بحرانی وجود دارد که فشار بخار سیال را تعیین می 
کند. این را می توان با استفاده از روش ]۳۹2 محاسبه کرد: 


(400.0) 51 . وم <<< 
9 466 


۳ ص دح 4 


نتیجه تا بیش از ۱۰ رقم دقیق است و با استفاده از برخی تکنیک‌های ریاضی که امکان محاسبه مستقیم فشار بخار را 
حل کننده نیوتن با دقتی که یک عدد ممیز شناور می‌تواند ارائه کند. پولیش می‌دهد: 


(400) ۲921 . 05ع/ (تاط1حطو 013 , 400) 921 1-09 <<< 
66-14 . 1 


چند روش دیگر مورد علاقه )۷1_98 و 826 ۷2 هستند که مایع اشباع و حجم مولی را محاسبه می کنند. ]152 که 
دمای اشباع را با یک فشار مشخص محاسبه می کند و 91_«ا که ضریب فوگاسیته اشباع را با توجه به یک دما 


(500) او 609۰۷ , (298۰15) 91 1 09.۷ <<< 
(۰0006827569 0 ,0۰0001303559) 

(0 ,۰ 101325) 1921 . وهوم <<< 

3 5 

(425۰0) 521 دطو. وهوم <<< 

0. 6 


۲۰-۱) کار با مخلوط ها 
استفاده از مخلوط 061100.605_101] ابتدا با استفاده از مخلوط هم‌مولار نیتروژن-متان در ۱۱۵ کلوین و ۱ مگایاسکال 


9 معادله پنگ-رابینسون نشان داده شده است: 


م [46.04۳5 ,945 . 33]عوعط , [190.6 ,126.1 ]عوم1 ,1۳۲6عظ ,۰0 ۲<2115) ۳872 < ووه <<< 
([[0۰0 ,0۰0289] , [0۰0289 ,0۰0] ]عو([3 رد.0 م5 .0و2 , [0.011 ,۰.04 0] 0۳6925 
605۰۷ ,1 ۵09.۷ <<< 

)3۰6587077706-05, 0۰00070676607( 

<<< 609.1۱093201 1]165 1, 6095 ۰, 1۱092011169 

)]838516۰99, 78350.27[, ]438108.61, 359993 ۰481[( 


تمام ویژگی های موجود در 00۳05 برای اشیاء 00۳05[۷116 نیز موجود است. اشیاء جدید 00۳0531 را می 


توان با متد 0] ایجاد کرد که کسرهای مولی جدید 75 و همچنین متغیرهای حالت جدید را می پذیرد. اگر یک ترکیب 
جدید 2٩‏ ارائه نشده باشد. ترکیب فعلی نیز برای شی جدید استفاده می شود. 


(5ه1عظ ,0 2<300) 0 .ووه <<< 

,0۰0111 ,۰04 0]عوووع۳ه , [0. 4604000 ,3394000.0]عوعط , [190.6 م,126۰.1]<وم۲) ۲۳12 
(100000۰0عظ ,۰0 300ع۲ ,0.5 مد ء 0و2 , [[0.0 ,0۰0289] , [0۰0289 م,0۰0] ] <و [ 11 
([9. م1.اعوه2 ,5ه1ع<ظ ,0 ۲<300) 609.10 <۶<< 


۳ ص دح هه 


,0۰0111 ,۰04 0] عحوووع۳ه , [01. 4604000 ,3394000.0]عوعط , [190.6 م,126۰.1]<وم۲) ۲۳۳12 
(100000۰0ع<ظ م۰0 300ع۲ ,[0۰9 م1 0و2 , [[0.0 ,0۰0289] , [0۰0289 م,0۰0] ] <و [ 11 
([6. م4.اعوه2 ,5ه1عظ ,)09.10 <<< 

,0۰0111 ,۰04 0] عحوووع۳ه , [۰.01 4604000 ,3394000.0]عوعط , [190.6 م,126۰.1]<وم۲) ۳۲۳۳12 
(<۷ ,۰0 100000عظ ,0.6 .0و2 , [[0.0 ,0۰0289] , [0۰0289 م,0۰0] ] <و [ 11 

([6. ,4 ءاعحوه2 ,1300.0 ,۰0 ۷<1) 605.1۲0 <<< 

,0۰0111 ,۰04 0] حوووع۳ه , [01. 4604000 ,3394000۰0]عوعط , [190.6 م,126۰.1]<وم۲) ۲۳126 
(0 ۰ ۷ ,300۰0ع۲ ,[0.6 م0۰4]عحوه2 , [[0.0 ,0۰0289] , [0۰0289 م,0۰0] ] <و [ 11 


ایجاد اشیاء جدید 01209۷11 با متد زیرمجموعه که فقط از برخی از اجزای اولیه مشخص شده استفاده می کند 


اسکاتفی تیگ 


م0۰00171] م [0.00061 م00 ,0۰0006] م [0.001711 ,0.00076 م,0.0]] < و11 <<< 
[ ۵۰.01 , 0۰00061 

م 0۰03311 ,0.1501 ,0۰8168 حوه2 , [540.3 ,50.4 ,469.7] <وم۲) ]۳۳ < ۳۳83 <<< 
,1۲2322۰29 , [2.136۳6 ,3.012۳6 ,369۳56 ,3] حوعظ , [0.349 ,0.305 ,249 0] <25 01969 
(روزاعو زنع ,۳<2101325:0 

([1,2]) ۲هوطداو . 283 <<< 

,۰305 0] <و ۳96۵02 , [۰.0 2736000 ,۰0 3012000] <وعظ , [40.3ظ ,۰.4 507] ع<وم1) ۳۲۳126 
,۰8193231441048 0]عو م2 ,0.01 ,0۰00061] م [0.00061 ,0۰0 ] حو 13 , [0.349 

(0 ,01325 1<ظ ,2322.29 , [ 8951 0.180676855 

([8. ,2. ]حور ,5ع1ع<ظ ,۲2۶500۰0 , [1,2]) 120961و , ۳۳83 <<< 

,۰.305 0] <19۵925ه , [2736000.01 ,۰.0 3012000]حوعط , [40.3 ,۰.4 907]<وم1۲) 1[1ظ 
۲2500۰0 م[0۰8 ,02]عحوه2 م[([0۰0 ,0۰00061] , [ 0.00061 ,0۰0 ] حو 11 , [0.349 
(0 .100000 <ظ 

([81. ,2.]<و2 ,[1,21]) 9۱10561 . ۳۲۳۵3 <<< 

,۰.305 0] <19۵95ه , [2736000.01 ,۰.0 3012000]حوعط , [40.3ه ,۰4 507]<وع1۲) 1[1ظ 
۲2۶322۰29 م[0.8 ,02]عوه2 م[([0۰0 ,0۰00061] , [0.00061 ,0۰0 ] حو 11 , [0.349 
(0 .01325 ۴-1 


همچنین می توان اشیاء 0005 خالص ایجاد کرد: 


<<< ۳۳83 , 012165 )( 

(۰0 01325 1عظ ,۶322.29 ,249. ۳98۵9220 ,۰0 3369000عوظ ,۰7 1۲6۵2469) 158 ] 

,۰3 540عع۲) 8 , (101325:0ع<ظ ,۲۶322۰:29 ,۰.305 ۳9۵9220 ,۰:0 3012000ععظ ,۰.4 ۲02507) ۳۲5٩‏ 
[ (0 01325۰ 1عظ ,29. ۲2۶322 ,349 .۵۳۵9220 ,۲022736000۰0 


۰/ ص وج ه 


دماء فشان تعداد مول 9 شتقات ۳3 مولی ضرایب فوگا بتی ورود به سب تم نیز با روش‌های 1 100۳019 


طل_دنطام‌صال عصل_عتطمصل و و2ل_فتطامرصل در دسترس هستند: 


(۲1۲) 01 و 010131 . 3 <<< 

]0 ۰029486952019, 0۰03514175794, 0۰0402818452753[ 

(۲1۲) 0۳ و 11و01 . 3 ح<< 

]-9.82537796-06, -9.81890930316-06, -9. 81225986-0 6[ 

(1 ۰2 ۴83) وصل و 1طمر01 . 3 <<< 

بت ,0۰0041532288] ,[0۰007300114797 ,0۰004153228837 ,0010590517 0۰-] ] 
,0۰02576802316 - ,0۰0073001147] , [0.0257680231- ,0۰016918292791 
[ 64621 0۰063291 

(1 ۰2 ۳83) 025 و 01۳0131 . 3 <<< 

] ]0۰0099380692, 0۰0151503498, 0.0182972357[ , ]-0.038517738, 0059589260, 
0.0684389901[ , ]-0.070571069, -0.103639207, -0.1411628301[ [ 


۴-۱) ویژگی های دیگر 


هش کردن 
با استفاده از متد ۲00061_50 می‌توان این دو شی را با یکدیگر مقایسه کرد و دید که آیا و1 پارامترهای مدل و 


اجزای یکسانی دارند یا خیر: 


,33۰94۳5 ] حوع‌ط ,[190.6 م126.1حوهم1 ,1۳6حط مط۲211) 126( < عوهن ظ۲ <<< 
( [[0.41,0] , [0,/0۰41] ]<و 1 م,[ظ.0 مدط.0]عحوه2 , [0.011 ,5<]0.04ع۳۵9ه , [04۳5 .46 


۰9485 33] حوع‌طظ , [190.6 م,109<2]126.1 ,1۳86حظ ,1۲2115) 9۳81126 < ووه له <<< 
([0[1 ,۰41 0] , [0,0۰41] ]<و لدع م[0.5 م,0.5] و2 , [0.011 ,0.04]حوووعوه , [ 46.045 


() 296.10006198 _کلظ۴و حح () 81و 10061. 6296 ظظ <<< 
6 ۳21 


با استفاده از متد 5]260_25 می توان مشاهده کرد که آیا حالت دقیق و مدل بین دو شی مختلف مطابقت دارند یا 


خیر: 


۳ ص و < 4 


,۰94۳5 33]حوع۴ , [190۰6 م,126.1]حوج1 ,6ظ۳لحط ,۲2116) ۲۳812 < 2ووه ۲۲8 <<< 
([[0 ,۰.41 0] م [0,0۰41] ]<و([1 م,[ظ.0 م5ء0] وه , [0.011 ,۰.04 0] وه , [ 46۰04۳5 


() 9 100061. 6962 ظظ حح () 91 10061. 0296 ۲8 <<< 
1۵ 
() 91816۵9 . 62962 _8ظ حح () 91 91216 .0296 ۲۲8 <<< 
6 ۳2۵1 


و در نهایت می توان با استفاده از روش _10881__ مشاهده کرد که آیا ده شی دقيقاً یکسان هستند» از جمله نتایج 
محاسبات کش شده: 


,)۰94۳۲5 33]حوعط , [190.6 ,9<]126.1ه1۲ ,6طظلاحط ,۲2115) 8126 < 03وعه ۲۳۲ <<< 
([[0 ,۰41 0] , [0,0۰41] ]<و لدع م[0.5 م0.5]عو2 , [0.011 ,0.04]حفووعهه , [ 46.045 


() 9162169 . 3ووه_ظظ حح () ۵91 91216 .0296 ۲۲8 <<< 
1 1 

(3ه‌ووه_5ظ) 1وقط حح (عوهه ۴8) 291ظ <<< 

1 1 

ورد 01 ح2طم1ع 0 . وووه 8 < <<< 

(3ه‌ووه_5ظ) 1وقط حح (عوهه 8ظ۴) 291ظ <<< 

۳21 6 


سریال سازی 
همه مدل‌های مکعبی ۳05 یک متد 25508 و یک 1701۳0501 را برای سریال‌سازی وضعیت شی برای انتقال از 
طریق شب؟ . ذخیره‌سازی روی دیسک و ارسال داده‌ها بین فرآیندها ارائه می‌کنند. 


موز 110۳0011 <<< 

م [11]حوم2 , [0.2975] <و9و1۵9ه , [30250001]عحوع‌ط , [501۰6]حوم1۲) ۳۶5۱2112 < ووه <<< 
([0۰460] حو120023 , [0.8634] <و 12۳0۳022 , [ ۰059104 0] عحو 12۳0021 ,16حظ ,۲-2299 
(() ط0وژ_ و . وهع) وصصتتل . طموژ < ۱0۶۶ _ط0وژ <<< 


((گناتاو_طموژ) 10208 , ط0و 3) طمو 10 , 60۳05112 < ومع ۵ <<< 
و << 05 ۳6۲۷ ]299611 <<< 


علاوه بر کتابخانه استاندارد 500[ می‌توان از کتابخانه‌های دیگر 502ز با طراحی استفاده کرد. ذخیره و بازآفرینی اشیا با 


کتابخانه وس .عنم پایتون نیز آزمایش شده است. این می تواند سریعتر از استفاده از 150 به مقدار باینری 


بودن داده باشد. 


۵-۱) تعادل مخلوط 


بر خلاف اجزای خالص.ء تعیین وضعیت تعادل برای مخلوط ها ساده نیست. الگوریتم های مختلفی مانند جایگزینی 
متوالی و کمینه سازی گیبس استفاده می شود. همه آن‌ها نیاز به حدس‌های اولیه دارند. که متعییو لا از مدل‌های 
ترمودیتامیکی ساده‌تر می‌آبند. در حالی که در عمل می توان تر کیب تعادل را برای یک مسئله لا فاز تعیین کرد در 


تتوری باید از یک الگوریتم بهینه سازی جهانی استفاده شود. 


جزئیات بیشتر در مورد این موضوع ۳ مین توانید در ماژول 161170.11250] بیابید 


۶-۱) استفاده از واحدها با معادلات حالت مکعبی 


برای استفاده از این اشیا نیز یک بسته بندی پینتی وجود دارد. 


۲ اط۳0د وا ذصتا . ۵110 مرگ << 
مک 1 [6 .190 م,126,.1]عحوم( ,1و ۶30*10 ,060 ۰0*1۰ ۲<۶400) 010۲ < و1720 <<< 
05]عو2 ,ووه 081 0116081 ۲9۵۰ [1 0.011 ,04 0] عوووموه ,معط ,1 [46./04۳5 ,۰.94۳5 33] عحوعظ 


(ووه‌لوه و01۳۵ , ۲۱۱ [ [0.0 ,0:0289] , [0۰0289 ,۰0 0] ] حو 11 


0.5 [ ۲۱ ۰ 01116۳۴531 0۴1695 


(۲1۲۷2۳199*) 8۳126 < واتصتا_ ووع <<< 
09 ,و (ع_ ۲ . وتاذصتا_ ووع <<< 
(<( ۲1۵1۳۷71۲ ,۰15 673) 11۷ 9029> ,<( 1۵161 7 ۲91 ,22۰.53858854 1:21111-۷7ه>) 


(2ظ ,۵6*1 1ع<ظ ,۰6 1300۰0۳1 16 < وووو <<< 


<<< 296۰.7 


>۵192۴111۷7 )0۰00249433879, ۱۲۳۵۲۵۲ ** 3 7 10۵161 < 


۳( ص دح ه 


۲ مقدمه ای بر مدل های ضریب فعالیت 

سیستم های تعادل بخارمایع و مایع-مایع می توانند رفتارهای متفاوتی داشته باشند. قانون رائول فقط می تواند 
وابستگی به دما و فشار را توصیف کند. بنابراین اغلب از یک ضریب تصحیح که وابستگی به ترکیب را اضافه می کند 
به نام "ضریب فعالیت " استفاده می شود. این یک رویکرد جداگانه برای استفاده از معادله حالت است. اما از آنجا که 
همبستگی فشار بخار مستقیم با ضرایب فعالیت استفاده می‌شود. می‌توان یک نتیجه با دقت بالاتر برای تعادل فاز به 
دست آورق: در حالی که این مدل‌ها اغلب «مدل‌های ضریب فعالیت» نامیده می‌شوند» در واقع یک پیش‌بینی برای 
انرژی اضافی گیبس هستند. ضرایب فعالیتی که برای تعادل فاز استفاده می شود از مشتق عدد مول جزئی انرژی 
گیبس اضافی بر اساس عبارت زیر بدست می آید: 


واع 0 ۶ 
01۳1 
٩ << 6۳ ]‏ 
۱ ) هه - ۱ 


۵ مدل ضریب فعالیت اساسی در ترمو وجود دارد: 


1( [۲ 

2( 

3( ۲ 

ممانا50 تقاناعع۴ (4 
۲0 (5 


هر یک از این مدل ها شی گرا هستند و از کلاس پایه 01005۳6655 که بسیاری از روش های رایج را ارائه می دهد. 
به ارث می برند. یک کلاس ساختگی دیگر که صفر اضافی انرژی گیبس و ضرایب فعالیت ۱ را پیش بینی می کند. به 
عنوان 1062150104100 موجود است. مدل انرژی اضافی گیبس معمولاً نسبتاً ساده است. تعدادی مشتق برای محاسبه 
سایر خصوصیات مانند ضریب فعالیت مورد نیاز است تا اين عبارات پیچیده تر از مدل واقعی به نظر برسند. در ادبیات؛ 
معمول است که یک مدل به طور مستقیم به شکل ضریب فعالیت بدون بحث در مورد مدل انرژی اضافی گیبس نشان 
داده شود. برای نشان دادن تفاوت. در اینجا بیان اثرژی گپیس مدل ,11و مدل ضریب فعالیت آن آمده است: 


11 ۹ 
84 ۳ " ۹ 
۱ 7 تا ۲ رت ۲ ‌ ۱ ۳ 7 زر ج) و1 1 ) ۷" 
تنس زز مزلم 1 ۱ 
مد وی سس ۷+ 3 - (:د)ا 


ص اب 


۲<1 


/ 1 2701 2 ۱ ۱ 2۳ ۰ عش ز 2۳0 َِ- 
ب ب ۳۳۳ سس 


۴/ ص دح 4 


مدل های ,۵1571 ۷۷1502 و 131017۸6 بیشترین استفاده را دارند. هر یک از آنها مبتنی بر رگرسیون هستند. هر 


یک از این مدل ها علاوه بر وابستگی به ترکیب. پارامترهای وابستگی به دما گسترده ای دارند. وابستگی‌های دمایی 
اعمال‌شده بایداجازه دهد که پارامترهای بسیاری از منابع دیگر در اینجا با آنها اتفاده شوند 


مدل 5682112750110108 مبتنی بر مفهوم پارامتر حلالیت است. با حجم مولی مایع و پارامترهای حلالیت یک مدل 
پیش بیتی است. وابستگی به دما را نشان نمی دهد. ضرایب ر گرسیون اضافی نیز می تواند با آن مدل استفاده شود. 
مدل »1/1۳۸ یک طرح مشارکت گروهی پیش‌بینی کننده است. در آن» هر مولکول به بخش های مختلف تقسیم 


می شود. این بخش ها دارای پارامترهای تعاملی با بخش های دیگر هستند. معمولا تکه تکه شدن با دست انجام نمی 
شود. یکی از ابزارهای آنلاین برای انجام اين کار 1001 احعصجعدعه ومنمتن عصنلمم 1۳۳85۲ است. 


۱-۳) ساختار شی 

شی 01005۳66595 چیزی در مورد تعادل فان فشار بخار يا روال فلاش نمی داند. دامنه آن به مقابله با انرژی اضافی 
گیبس محدود است. به دلیل این مدولار بودن» یک شی 01005۳70655 مقداردهی اولیه طراحی شده است تا در یک 
آرگومان به معاذلات مکمب.حالضی که از ارتی گیسن اخافی ماندد ۹3۴ اسفاده می کند» ارسال شوذ. 


مکان دیگری که این اشیاء مورد استفاده قرار می گیرند در اشیاء 0010 ]01009106655 است که قطعات را به هم 
نزدیک می کند تا یک فاز ترمودینامیکی (عمدتا) سازگار ایجاد کند که الگوریتم های فلاش می توانند با آن کار کنند. 
این ماژولار بودن به مدل‌های اضافی گیبس اجازه می‌دهد تا در هر مکانی نوشته و استفاده شوند. بنابراین مدل ۳61 
با خوشحالی به یک شی 1/1۳۸ اجازه می‌دهد که مانند ۷۲۳1۶ پیکربندی شده باشد. 


۲-۲ مثال 11۲۸ 
فول ها بکه‌سشالن یی شیب متا کت گروفی است. گضنا اسضاه آ کیتوی 6 هر مولکول یه 


تعدادی «گروه» و «شمارش» آن‌ها کار می‌کند. 


موه مک کلات: فده را با ادها ۱۳۳۵ مر کردم اس یک سیستی فوتایی ساخهاه سل 


2 در ۲ با «هگزان و بوتانون-۲ در زیر نشان داده شده است: 


0 ,۱7۳96 مط۳1] اه1۳ ۶20 ده ۰ ۵110 مرگ << 

م [ 18:1 م2:1 ,1:1 م,(2:4 ,1:2 ] حوصدا۵۳۲9۲ظه) وصتا۲0وطناو 10۳8 7۵]۳20] < لت <<< 
,02121۳1 16120۲10 , 0احطم 7۵91 ر[5. ,5 0]حوعد رک 1<60+273 
100۴0120920556 و 


۱۵ ۱ ص غ< ه 


راه حل ارائه شده توسط ۳135۲ دارای مقادیر ضریب فعالیت [۰۱,۴۲۸ ۱,۳۶۵] است که با موارد محاسبه شده توسط 
شی 1/(1۳۸0 مطابقت دارد: 


<<< 0۳ . 0211۳125 )( 
]1 ۰. 4276025835, 1۰3646545010 [ 


بسیاری از ویژگی های دیگر نیز پیاده سازی شده اند که تعدادی از آنها در زیر نشان داده شده است: 


(0 0۳2 6.0288 , () 8۳_ عم . هم ر () و .هع <<< 
(0:00380070204- ,0۰206721488 ,923.641197) 

() ۳_ 085 . ظ6 , () ۳_طتله . قم , () ظ.ق , () ظ1.ت <<< 
(0:0038007020460 ,1:266203886 ,0۰2067214889- ,۰77193363 854) 


توجه داشته باشید که متغیرهای ۳1۳ و 0۳50 حاوی پارامترهای تعامل واقعی هستند. هیچ کدام با کلاس هاردکد 
نشده آند. بنابراین کلاس می تواند برای رگرسیون استفاده شود. پارامتر نسخه کنترل می کند که کدام نوع از »۱01۳۸ 
استفاده شود زیرا تعداد کمی از آنها وجود دارد. مدل‌های مختلف 171۳۸6 پیاده‌سازی شده شامل »111۳۸ اصلی. 
0 دور تموند. کل ۷۲۳ مهوضا/0عصس و ]16 11۳۸ است. پارامترهای تعامل برای همه مدل 
ها نیز گنجانده شده است. اما آرگومان نسخه به فایل های داده متصل نیست. 


برای راحتی. تعدادی از تکه تکه شدن مولکول ها با کد 1۳۸( توزیع می شوند. تمام قطعه بندی ها از طریق 
پورتال آنلاین 0۳857 که در آن فایل های ساختار مولکولی را می توان ارسال کرد. به دست می‌آید و این مزیت را 
دارد که آنجه اراته می شود بدون ابهام است. هیچ نگرانی در مورد اعداد ۸5) وجود ندارد» مانند اينکه چگونه گرافیت 
و الماس شماره )2۸٩‏ متفاوتی دارند در حالی که یک عنصر هستند يا هوا با وجود مخلوط بودن. دارای شماره 6۸۵5 
است. بر این اساس» فهرست موجود در این فایل‌های داده توزيع‌شده. کلیدهای 1001 هستند که می‌توان آنها را از 
وعزکنامعل1. علمءتصعطه پا در مکان‌های مختلف آنلاین دریافت کرد. 


<<< 110۴001 21۲0 . 

() ۲1355 2553198106۳185 ت01 1 1080 .17۸31۴20 ۰ 161100 <<< 
(695ع1طو 2553 17۱11۳2۸0 ۲فطنار , ۴26 231 ۰ 1261100) طه 1 <<< 
29946 
09۳06۳95 2951 17۵11۳۸0 (1۳1 011۳ ۲ وظنار , ۴26 ۰7231 1261100) طه 1 <<< 
29971 
(5اط6عتطو 51و تطوظ ‏ ۲ فوطنان , ۴26 ها , 1261100) طه 1 <<< 
320034 

1 هه _ اجه هرود وه تفه مر <<< 


۷ 180۳81 . ( ۲۲۵1۱0۵86۲ ) 691 ت۵طه _طمزوهو << 
5۵-۲ ۲۲۸۵۷ ۱۲۲ -۵2)۳۴۹۵(۲۲۱۲ ۱۷6۲۱/۷2۵ ۲ 


۳ ص دح ه 


-۵250(۲۲۲ ۲۲۷2۵۲۷/۷۸ ] ۲5 منطو 1 ووه 17۸۲1۳۸0 1۳۱ 101۲1۳ فرظ , ۴26 صنا . ۵1۲80 <<< 
[ ۲ 1111۳۳۳۸۵۷۹۵1۲ 
ال بل میت ع 1 


لطفا توجه داشته باشید که عدد صحیح شناسایی در این عناصر (اصنمی :نامع ۳ در نسخه‌های مختلف 
0 بیکسان نیست. 


۲-۳) نکاتی در مورد عملکرد 

مقداردهی اولیه شی برای اولین بار یک عملیات با کارایی بالا نیست زبرا باید بررسی های خاصی انجام شود و 
ساختارهای داده تنظیم شود. برخی از قطعات معادلات مدل اضافی گیبس ممکن است فقط به دما يا ترکیب بستگی 
داشته باشند. نه اينکه به هر دو بستگی داشته باشند. هر مدل روش 01_55 را پیاده سازی می کند که باید برای 
ایجاد یک شی جدید در دما و/یا ترکیب جدید استفاده شود. طراحی شی برای محاسبه خواص و غیرقابل تغییر بودن 
است: محاسبات در دماها و ترکیبات جدید در یک شی جدید انجام می شود. همچنین توجه داشته باشید که رشته 
_ 160۲ __ برای هر مدل طوری طراحی شده است که امکان بازسازی بدون اتلاف مدل را فراهم می کند. این در 
هنگام ساختن موارد آزمایشی بسیار مفید است. 


([9. ,1 ]عفد ,۲-400.0) ۶۵ ۲ 2.20 <<< 
,3.856]عوو , [3.2479 ,4998000000000005 4] 29 , [0۰9 ,0.1] 9 ,۰0 ۲-400) 0111۳۸6 
ر 11 ۱0 11 ۱۸ 11 مما] ]عفن ر[188 0۵54 948: 0 ]وه [2:876 

0۰01 ,26.76 ,26:76] , [476:4 ,0:0 ,0.0] [4: 6 ,0۵ ,0:0]])عه‌طه_ زوم 
و0 وه ای ]۱ 0 م0 0 0 مه یل ره لا ۲0:0۶ )۲ 
(ون وه م( 00 منم و0 و00 1:0 0 :0] 


هنگام کار با تعداد کمی از مولفه ها (۵ يا کمتر). ۳۷۳۷ بهترین عملکرد را ارائه می دهد و استفاده از مدل با لیست 
های ۳:07 به عنوان ورودی سریع ترین راه برای انجام محاسبات حتی در ۳۷)00) است. اگر با بسیاری از مولفه 
ها کار کنید یا اگر آرایه های :۱۷۵ به عنوان ورودی و خروجی مورد نظر هستند. آرایه های ۱۷:0 می توانند به 
عنوان ورودی ارائه شوند. این تأثیر منفی بر عملکرد خواهد داشت مگر اینکه از رابط 007002 استفاده شود: 


صص وج تا 110۳00 << 

تاد . معط اطم100 <<< 

<<< 1 - 3 

<< 71 - 25.0 + 2315 

([0.1818 ,0.0909 ,0۰7273]) 1۳,۵12 < ویر <<< 
([3.1878 ,2.1055 ,92.]) 102112۷7 < 25 <<< 
([2.4 ,1.972 ,1.4]) 0۰2112۷ < ون <<< 


۱/۷ ص دح ه 


طد ند طو ۲( [0۰0]]) 1۷ج 9ظ << واه < 21895 < نوتاه < 12۵1790 < 12192 <<< 

12۳06 )5( [( 

م1325.1-] م[91.532 ,لا ,۰06 318] , [309.64- ,526۰.02 م,0]]) ۲۵۰۵1۳2۷ < طویاه1]2 <<< 
را 7ص فص ل وت 

(۲طوناه] ,واه رناوناه] ,نوتاه ,طفتاه] مطوتاهع]) < ۲سانآنظ2۳ <<< 

0 1۳00۴ ۵0 ورگ << 

(دطنآن‌طحع نانک , حول , وطحوط ,2:95 , 11) 19۵170 < 6۳2 <<< 

<<< 6۳2 . 02۳۳125 )( 

2۳11۳12۷7 )] 1.5/7039333, 0۰29482416, 18 .11432905[ ( 


رابط 1001002 سرعت کد را افزايش می دهد و امکان محاسبات با ده ها مولفه را فراهم می کند. رابط 000102 نیاز دارد 
که همه ورودی‌ها آرایه‌های 001001 باشند و همه خروجی‌های آن نیز آرایه‌های 0101002 هستند. 
( نان طحع نانک , وحوه ,و حوط ر 5(حوعد , 11) لام 1 لا . و9متاظ . ۵1100( << سای <<< 


<<< 6۳3 . 27۳125 )( 
21۳12۷ )] 1.5/039333, 0.29482416, 18 ,11432905[( 


به عنوان نمونه ای از مزایای عملکرد. محاسبه گامای »۳۸ ۲۰۰ جزء. ۱۰,۶ میلی ثانیه در 0۳۷00۴ و ۳۲۱۸ 
میکروثانیه زمانی که توسط 1۷01002 شتاب می گیرد. طول می کشد. در این مورد ۳۳۷ ۶۶۴ میکرو ثانیه طول می 
کشد. هنگامی که همان معیار با ۱۰ مولفه انجام می شود. محاسبه در 0۳۷008 ۲۸۷میکروثانیه. با 20۳002 
۶میکرو ثانیه و با ۲۳ ۳۶,۲ میکرو انیه طول می کشد. 


استفاده مجدد از روش 01_55 می تواند بسیار مهم باشد. برای »۰11۳۸ چندین اصطلاح فقط به دما بستگی 
دارد. اگر محاسبه ۲۰۰ مولفه با مواردی که قبلاً محاسبه شده تکرار شود زمان‌بندی‌ها در 0۳02 ۲۶.۳ میکرو 
ثانیه» با 20۳002 ۱۲۷میکرو ثانیه و با ۳۷۳۷ ۲۵ ۱میکروثانیه است. 


۲-۴) ویژگی های دیکر 


ضرایب فعالیت محدود کننده رقت نامحدود را می توان با فراخوانی به طمناباانه_عتتصتلص_عمصصعع به دست آورد. 


<<< 0۳ . 211025 18۴181 6 0311063 05 )( 
]3 ۰5659995166, 4۰.32849696[ 


۱/۸ ص غ< ه 


همه مدل‌های ضریب فعالیت یک متد 29507 و یک 1۳01050۳0 را برای سریال‌سازی وضعیت شی برای انتقال از 
طریق شبکه. ذخیره‌سازی روی دیسک و ارسال داده‌ها بین فرآیندها ارائه می‌کنند. 


0 1231 501 10621 امد معطاا مر << 

موز 110۳001 <<< 

([4. موه م2 م1 عوع< ,۲<300۰0) 10۵215011110۳1 << 10061 <<< 

() 501 3 و , ۲۵061 < ۷1۵7 ط0ود <<< 

(۲۷16۷_طهو ژ) دول . طموژ < او طووژ << 

((او_ طهو 3) 10208 , هو ژ) طمو 101 , 0 1013 10621501 < 00۵۷ 0061 <<< 


<<< 299611 10061 000۷ << 1 


علاوه بر کتابخانه استاندارد 500[ می‌توان از کتابخانه‌های دیگر 500 با طراحی استفاده کرد. ذخیره و بازآفرینی اشیا با 
کتابخانه دوس .عنم پایتون نیز آزمایش شده است. این می تواند سریعتر از استفاده از 150 به مقدار باینری 
بودن داده باشد. همه مدل‌ها دارای یک روش _ 250 __ هستند که می‌توان از آن برای مقایسه مدل‌های مختلف 
ستفاده کرد تا ببينيم آیا کاملاً یکسان هستند (از جمله اینکه چه مقادیری قبلاً محاسبه شده‌اند). 

آنها همچنین دارای یک متد 1710061_025۳ هستند که می‌توان از آن برای مقایسه مدل‌های مختلف استفاده کرد تا 
ببینیم آیا پارامترهای مدل یکسانی دارند یا خیر. 

آنها همچنین دارای یک متد 216_250) هستند که می‌توان از آن برای مقایسه مدل‌های مختلف استفاده کرد تا 


ینمی آیا دماء ترکیب و پارامترهای مدل یکسانی دارند یا خیر. 


۵-۲) هوبت های ضریب فعالیت 


مجموعه ای از معادلات مفید به شرح زیر است. برای اطلاعات بیشتر به منابع مراجعه شود. هیچ منبعی حاوی همه 


۳ ص و < 4 


منابع 


۰و 20 وموعت) ۵۶ وع۳۲۵0۵۵۲ ۲۳ باامصجمن 0۵۰ .۲ ول 4صع رماتصعننه:۳ ,۱ صط10 ر.ط ومتاظ رعصتانظ [1] 


0۰ ۳۲0۶65۹101861 ۱۷]۵0۲۵0۷-۳۲111 ۷۵۲۱۰ ۱۱۵۷ ,وملات0ع طاد 


تمه رل مصط تمط ز. آهمتمصمطی) ۵۲۵۲۰ ممععتا [ مره ,امک [مماقظ رتمک آممطم۱۷]1 بممع‌عنا ٩‏ معصتلطمصرت [2] 


۰ ,088 عک ۷۷۱۱۵ ام مصمتاها نگ وقععمعظ ۲۵۲ 


,۷۵۷۰ 6010100۰ 2۵0 .۵۲5م1مدظ لدمتصمصن ۲۵۲ وتات وظ همه مره لهمتووطظ ع0 آع0 روته ۵۷ [3] 


2012 


تون مصمتات0ه 2۳00 .ومنصصههل0مصصمط 1 مصیممصلمص تمهت ۲6۲۵006۲۵۲۷ هن آنهن 200 و.ل متامتااظ [4] 


520016 ۳۲۷6, ۲: ۲۲۵۵۲۱۵6 ۳۲۱۵1۱۷ ۰ 


5۰ مصمصصم ماع ]6-۲ ع0ع با بعمتممصتم مط لهمتصصمطی م1 ۳0111012 مفقطاظ ,۱۷۲ مها ۲۳ رفهله ۷۷ [5] 


۳۰( صرح 4 


